Global phylogeography of ridley sea turtles (Lepidochelys spp.): evolution, demography, connectivity, and conservation.

Globally distributed marine taxa are well suited for investigations of biogeographic impacts on genetic diversity, connectivity, and population demography. The sea turtle genus Lepidochelys includes the wide-ranging and abundant olive ridley (L. olivacea), and the geographically restricted and 'Critically Endangered' Kemp's ridley (L. kempii). To investigate their historical biogeography, we analyzed a large dataset of mitochondrial DNA (mtDNA) sequences from olive (n = 943) and Kemp's (n = 287) ridleys, and genotyped 15 nuclear microsatellite loci in a global sample of olive ridleys (n = 285). We found that the ridley species split ~ 7.5 million years ago, before the Panama Isthmus closure. The most ancient mitochondrial olive ridley lineage, located in the Indian Ocean, was dated to ~ 2.2 Mya. Both mitochondrial and nuclear markers revealed significant structure for olive ridleys between Atlantic (ATL), East Pacific (EP), and Indo-West Pacific (IWP) areas. However, the divergence of mtDNA clades was very recent (< 1 Mya) with low within- clade diversity, supporting a recurrent extinction-recolonization model for these ocean regions. All data showed that ATL and IWP groups were more closely related than those in the EP, with mtDNA data supporting recent recolonization of the ATL from the IWP. Individual olive ridley dispersal between the ATL, EP, and IN/IWP could be interpreted as more male- than female-biased, and genetic diversity was lowest in the Atlantic Ocean. All populations showed signs of recent expansion, and estimated time frames were concordant with their recent colonization history. Investigating species abundance and distribution changes over time is central to evolutionary biology, and this study provides a historical biogeographic context for marine vertebrate conservation and management. Supplementary Information: The online version contains supplementary material available at 10.1007/s10592-022-01465-3.

Age-specific growth and maturity estimates for the flatback sea turtle (Natator depressus) by skeletochronology.

To address a major knowledge gap for flatback sea turtles (Natator depressus), a species endemic to Australia and considered 'Data Deficient' for IUCN Red List assessment, we present the first-ever skeletochronology-derived age and growth rate estimates for this species. Using a rare collection of bone samples gathered from across northern Australia, we applied skeletochronology and characterized the length-at-age relationship, established baseline growth rates from the hatchling to adult life stages, and produced empirical estimates of age-at- and size-at-sexual-maturation (ASM, SSM). We analyzed humeri from 74 flatback sea turtles ranging in body size from 6.0-96.0 cm curved carapace length (CCL), and recovered from Western Australia (n = 48), Eastern Australia (n = 13), central Australia (n = 8; Northern Territory n = 3, the Gulf of Carpentaria n = 5), and unknown locations (n = 5). We identified the onset of sexual maturity for 29 turtles, based on rapprochement growth patterns in the bones. Estimates for ASM ranged from 12.0 to 23.0 years (mean: 16.3 ± 0.53 SE), SSM ranged from 76.1 to 94.0 cm CCL (mean: 84.9 ± 0.90 SE), and maximum observed reproductive longevity was 31 years for a 45-year old male flatback. Growth was modeled as a smoothing spline fit to the size-at-age relationship and at the mean SSM (84.9 cm CCL) corresponded with a spline-predicted maturity age of 18 years (95% CI: 16 to 24), while mean nesting sizes reported in the literature (86.4 to 94 cm CCL) corresponded to estimated ages of 24+ years. A bootstrapped von Bertalanffy growth model was also applied and showed consistencies with the spline curve, yielding an estimated upper size limit, Linf, at 89.2 ± 0.04 cm (95% CI: 85.5 to 95.9 cm) with the intrinsic growth rate parameter, k, at 0.185 ± 0.0004 (0.16 to 0.22); at the same mean SSM (84.9 cm CCL) the estimated ASM was 16.3 ± 0.05 years (95% CI: 12.8 to 27.7 years). Lastly, four of the samples analyzed were collected from deceased adult females that had previous sizes known from on-going mark/recapture studies at nesting sites in Western Australia. The paired CCL data (measured at nesting and back-calculated) did not significantly differ (p = 0.875). This first skeletochronology study for flatback sea turtles generates valuable empirical estimates for ongoing conservation and management efforts.

Geography best explains global patterns of genetic diversity and postglacial co-expansion in marine turtles.

For many species, climate oscillations drove cycles of population contraction during cool glacial periods followed by expansion during interglacials. Some groups, however, show evidence of uniform and synchronous expansion, while others display differences in the timing and extent of demographic change. We compared demographic histories inferred from genetic data across marine turtle species to identify responses to postglacial warming shared across taxa and to examine drivers of past demographic change at the global scale. Using coalescent simulations and approximate Bayesian computation (ABC), we estimated demographic parameters, including the likelihood of past population expansion, from a mitochondrial data set encompassing 23 previously identified lineages from all seven marine turtle species. For lineages with a high posterior probability of expansion, we conducted a hierarchical ABC analysis to estimate the proportion of lineages expanding synchronously and the timing of synchronous expansion. We used Bayesian model averaging to identify variables associated with expansion and genetic diversity. Approximately 60% of extant marine turtle lineages showed evidence of expansion, with the rest mainly exhibiting patterns of genetic diversity most consistent with population stability. For lineages showing expansion, there was a strong signal of synchronous expansion after the Last Glacial Maximum. Expansion and genetic diversity were best explained by ocean basin and the degree of endemism for a given lineage. Geographic differences in sensitivity to climate change have implications for prioritizing conservation actions in marine turtles as well as for identifying areas of past demographic stability and potential resilience to future climate change for broadly distributed taxa.

Phylogeography, Genetic Diversity, and Management Units of Hawksbill Turtles in the Indo-Pacific.

Hawksbill turtle (Eretmochelys imbricata) populations have experienced global decline because of a history of intense commercial exploitation for shell and stuffed taxidermied whole animals, and harvest for eggs and meat. Improved understanding of genetic diversity and phylogeography is needed to aid conservation. In this study, we analyzed the most geographically comprehensive sample of hawksbill turtles from the Indo-Pacific Ocean, sequencing 766 bp of the mitochondrial control region from 13 locations (plus Aldabra, n = 4) spanning over 13500 km. Our analysis of 492 samples revealed 52 haplotypes distributed in 5 divergent clades. Diversification times differed between the Indo-Pacific and Atlantic lineages and appear to be related to the sea-level changes that occurred during the Last Glacial Maximum. We found signals of demographic expansion only for turtles from the Persian Gulf region, which can be tied to a more recent colonization event. Our analyses revealed evidence of transoceanic migration, including connections between feeding grounds from the Atlantic Ocean and Indo-Pacific rookeries. Hawksbill turtles appear to have a complex pattern of phylogeography, showing a weak isolation by distance and evidence of multiple colonization events. Our novel dataset will allow mixed-stock analyses of hawksbill turtle feeding grounds in the Indo-Pacific by providing baseline data needed for conservation efforts in the region. Eight management units are proposed in our study for the Indo-Pacific region that can be incorporated in conservation plans of this critically endangered species.

Determining microsatellite genotyping reliability and mutation detection ability: an approach using small-pool PCR from sperm DNA.

Microsatellite genotyping from trace DNA is now common in fields as diverse as medicine, forensics and wildlife genetics. Conversely, small-pool PCR (SP-PCR) has been used to investigate microsatellite mutation mechanisms in human DNA, but has had only limited application to non-human species. Trace DNA and SP-PCR studies share many challenges, including problems associated with allelic drop-out, false alleles and other PCR artefacts, and the need to reliably identify genuine alleles and/or mutations. We provide a framework for the validation of such studies without a multiple tube approach and demonstrate the utility of that approach with an analysis of microsatellite mutations in the tammar wallaby (Macropus eugenii). Specifically, we amplified three autosomal microsatellites from somatic DNA to characterise efficiency and reliability of PCR from low-template DNA. Reconstruction experiments determined our ability to discriminate mutations from parental alleles. We then developed rules to guide data interpretation. We estimated mutation rates in sperm DNA to range from 1.5 × 10(-2) to 2.2 × 10(-3) mutations per locus per generation. Large multi-step mutations were observed, providing evidence for complex mutation processes at microsatellites and potentially violating key assumptions in the stepwise mutation model. Our data demonstrate the necessity of actively searching for large mutation events when investigating microsatellite evolution and highlight the need for a thorough understanding of microsatellite amplification characteristics before embarking on SP-PCR or trace DNA studies.

DNA barcodes for globally threatened marine turtles: a registry approach to documenting biodiversity.

DNA barcoding is a global initiative that provides a standardized and efficient tool to catalogue and inventory biodiversity, with significant conservation applications. Despite progress across taxonomic realms, globally threatened marine turtles remain underrepresented in this effort. To obtain DNA barcodes of marine turtles, we sequenced a segment of the cytochrome c oxidase subunit I (COI) gene from all seven species in the Atlantic and Pacific Ocean basins (815 bp; n = 80). To further investigate intraspecific variation, we sequenced green turtles (Chelonia mydas) from nine additional Atlantic/Mediterranean nesting areas (n = 164) and from the Eastern Pacific (n = 5). We established character-based DNA barcodes for each species using unique combinations of character states at 76 nucleotide positions. We found that no haplotypes were shared among species and the mean of interspecific variation ranged from 1.68% to 13.0%, and the mean of intraspecific variability was relatively low (0-0.90%). The Eastern Pacific green turtle sequence was identical to an Australian haplotype, suggesting that this marker is not appropriate for identifying these phenotypically distinguishable populations. Analysis of COI revealed a north-south gradient in green turtles of Western Atlantic/Mediterranean nesting areas, supporting a hypothesis of recent dispersal from near equatorial glacial refugia. DNA barcoding of marine turtles is a powerful tool for species identification and wildlife forensics, which also provides complementary data for conservation genetic research.

The genetic structure of Australasian green turtles (Chelonia mydas): exploring the geographical scale of genetic exchange.

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Phi = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.